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Abstract 

We study cosmological aspects of braneworld models with a warped dimension and 
an arbitrary number of compact dimensions. With a stabilized radion, a number of 
different cosmological bulk solutions are found in a general case. Both one and two 
brane models are considered. The Friedmann equation is calculated in each case. 
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1 Introduction 



The study of extra dimensions has [Q , in addition to the novel features of particle phenomenology, 
inspired much research on the cosmological aspects of extra dimensional models. Especially the 
study of brane cosmology has recently been an active field [2]. Already very soon after the pioneering 
papers on brane world models the peculiar cosmological aspect of these constructions were 
realized Later on, a number of papers have studied the cosmological bulk solutions of RS-type 
models, |£|-[Z|. 

An important aspect of any cosmological brane world model is the question of radion stabiliza- 
tion, which in the RS-model corresponds to the stabilization of the distance between the two three 
branes. At the same time one wishes to solve the hierarchy problem introducing a fundamental, 
D -dimensional gravity scale M* which is related suitably to the Planck mass Mpi, at least at the 
present time. Most interesting would be if the fundamental scale M* is near 1 TeV and hence in 
reach of the future colliders. 

In addition to the many studies on ADD [S] and RS [3] -models, a number of papers have 
combined the two so that in addition to the warped dimension, a compact extra dimension is 
also present, see e.g. |§]-|13|. The stabilization of six dimensional models have been discussed in 
|14| I15j. Combining the two models can be seen as a quite natural scenario since in any case we 
expect to have a number compact dimensions coming from string theory. Obviously these compact 
dimensions do not have to be 'large' but such a possibility is worth studying. Thus of the D-l spatial 
dimensions three are the usual observed dimensions, one is a Randall-Sundrum like dimension and 
D-5 are compact, but unconstrained dimensions. Hence the dimension of the brane is D-2. Because 
the Standard Model particles are free to propagate in the extra compact dimensions, their size are 
constrained by experiments. In particular in our models, which have the extra dimension topology 
of J x (S" 1 )^- 5 ), where I C IR, the radii of S 1 factors are constrained to be smaller than about 1 
TeV" 1 . 

In this paper we are interested in the cosmological properties of brane world models, where 
the number of the extra dimensions is two or more. It is organized as follows: First we study 
the D-dimensional (D > 5) Einstein's equations in order to find some cosmological bulk solutions. 
Writing the D-dimensional jump conditions we solve the Friedmann equation in each case. Then 
we concentrate on the 6-dimensional models and their properties and finally present conclusions. 

2 Einstein's equations in D-dimensions 

In order to write the Einstein's equations in D-dimension, we first define the coordinates 1 as: 

x M = (t,x 1 ,x 2 ,x 3 ,z,0i, ...,0d_ 5 ). (1) 

Due to the symmetries, the metric tensor can always be chosen to have diagonal form with bi = 
bi(t, z). Hence the line element can be written as 

D-5 

ds 2 = rj(t, zfdt 2 - R(t, z) 2 dx i dx i - a(t, zfdz 2 - ^ b t (t, zfdd 2 . (2) 

i=i 

With this form of the metric tensor, the Einstein's equations can be written compactly as 

Gab = -k 2 Tab, (3) 

: In this paper we follow a convention where capital roman letters denote the full D-dimensional space- 
time, small roman letters the three dimensional space and small greek letters the four dimensional space-time, 
e.g. A = 0, D-l, a= 1,2,3, a = 0,...,3. 
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where k is the D-dimensional gravitational constant, Gab is the Einstein's tensor in D-dimensions 
and Tab is the D-dimensional energy-momentum tensor. It includes both the bulk and the brane 
energy-momentum tensors. In the bulk the energy- momentum tensor is assumed to have the form 

Tf = dmg(-A B ,-A B ,-A B ,-A B ,T*(t,z),...,TE:l(t,z)). (4) 

In other words, we have a normal cosmological constant in the bulk but the cosmological constants 
in the compact dimensions are allowed to be functions of time and the ^-coordinate. As it was 
already seen in jE], this is necessary in order to satisfy the continuity (and Einstein's) equations in 
the bulk. 

We now look for solutions with a stable radion field, i.e. a = 0. We are free to choose it to have 
a constant value in z-direction by redefinition of z-coordinate, thus simply a(t, z) = 1. We assume 
further that the compact dimensions are also stabilized, bi(t,z) = bi(z). Non-zero components of 
Einstein's tensor are now: 

D ~ 5 b'b' D ~ 5 b" D ~ 5 b'- R' R' 2 R" R 2 

G °° = .^ib J + ^i +3 ^b l R +3 jp +3 ii- 3 ^& 

D-5 yiyi D-5 yi D-5 u / D-5 u p / / p / p /2 

^hbj {^h ^b iV ^hR rjR R 2 

ri" R" nR R 2 R „ 

- + 2— + 24- --_ -2_ (1 = 1,2,3 
rj R j] 6 R rj z R z r] z R 

~ i£>i b i b j vh b * R h h t \ R R 2 v 3R v 2 R 2 v 2 R [) 

•°~ 5 h'h' D ~ 5 W P' D ~ 5 U ->V P' Q P' 2 m" 

R" T] R R 2 R . 

3 ll +3 ^R- 3 ^ " 3 ^fR {n = 5 ' 

rj'R R' 

G 4 — -3 — ^"+3— 
7] R R 

As we see, the Einstein's tensor does have a non-diagonal component G04 whereas the energy- 
momentum tensor is diagonal. Thus as the first step we solve The equation G04 = by writing 
r/(t, z) in the form 

rj(t,z) = X(t)R(t,z), (6) 

where A is a arbitrary positive time dependent function. This is naturally the same solution that 
was found in the 5-dimensional case [HI El • 

Using the solution © for ri(t, z), the non-zero components of become rewritten as 

o h'h' D ~ 5 h" D ~ 5 h' P' P' P" 

l^J>\ 1=1 1=1 

, 1 b% D ^b'l f R> R'^Vi R' 2 

X R ^ . b ' h ' R R h b * R2 

p" i, p' p' p" 

it a it it it ,, . ,_. 

A 2 i? 2 i^bibj R R ~[h R 2 \3RR RR 
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o D-5 UlJ D-5 in r>i p/ D— 5 »/ tv2 „// 

A R 1 R' R" , _ , 

3 r + 3 ^ + ^- (n = 5,...,D-l). 

A 3 i?i? RR R K 1 

The Einstein's equations which are directly related to the 4-dimensional space-time properties 
via couplings to 4-dimensional matter are the equations for Gq and G\. In contrast the other 
components of the Einstein's tensor are related to the unknown components T£, ...,t5Z\ of the 
energy-momentum tensor and remain, at this stage, less constrained. 

Besides the Einstein's equations, the equations to be satisfied include also the continuity equa- 
tions in D-dimensional space-time, T AB B = 0. However, in the bulk only one of these continuity 
equations is non-trivial giving 

T AB B = Tf + (T 4 4 - A B )(3^ + 5) + E( T 4 4 ~n) h f = 0- (8) 

We remind the reader, that the continuity equations are automatically satisfied due to Bianchi 
identities whenever the Einstein's equations have been solved. 



3 Cosmological bulk solutions 

To reach our goal to study the cosmological properties of the brane models we have to solve first 
the Einstein equations © in the bulk. The whole set of equations nor even the first four are hardly 
solvable in general. However, our special choices allows us to relate these to each other. 

Thus in order to solve the Einstein equations we simplify the set of equations (J2J) further. We 
first note that since r)(t, z) is given by ((Sj). we can write G\ in terms of G\: 

Gi = G° + ljd t G° (no sum), (9) 

hence in the bulk, as long as is constant, we only need to look for solutions of the first Einstein's 
equation Gq = k 2 Ab in order to solve the first four Einstein's equations. 

Even now, the general solution of the equation Gq = k 2 Ab is not known. We can, however, try to 
solve it in two simple special cases: We can use either exponential Ansatze where it is assumed that 

2 

bi(z) = exp(/cjz) or a cosh -Ansatz, which we call the bowl-model, where bi{z) = 6j(0) cosh d +! (kz) 
and d = D — 5 is the number of compact extra dimensions. These two choices carry some essential 
and interesting properties as will be seen. We first study the exponential model and then the 
bowl-model. 



3.1 Exponential models 

We first assume the exponential Ansatze for the scale factors of the compact dimensions, bi(z) = 
fei(0) exp(/cjz). The (OO)-component of the Einstein's equations in the bulk is then 

o D-5 jji D—5 r>r o// 

" + E hkj + 3^ J> + 3(^) 2 + 3^ = . 2 A B , (10) 

i>j=l i=l 

which is easily solvable. The nature of the solution depends, however, on the relations between the 
parameters of the equation. Indeed, there are four distinct classes of solutions, in addition to the 
trivial case k{ = Vi, to equation Eq. (fTU|). depending on the values of the constants hi. 
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By denning the combinations of the constants fcj as A = J2i>j=i 

kikj and B = £f ~ 5 h the 

solutions to Eq. (fTU|) are given by 

A) B 2 ^^(A-k 2 A b )^0: R 2 (t,z) = 3/(X 2 (t)(A - k 2 A b )) + 

di(i) exp (|(-3B + ^9B 2 - 24(4 - k 2 A b ))) + 

d 2 (i) exp (|(-3S - ^952 - 24(4 - K 2 A B )) 

B) A = k 2 A b , B / : R 2 (t,z) = 2z/(X 2 (t)B) + d^t) exp(-Bz) + d 2 {t) (11) 

C) A = k 2 A b , B = 0: R 2 (t,z) = di(t) + d 2 (t)z + z 2 /X 2 (t) 

D) B 2 = ^(A-k 2 A b )^0: R 2 (t,z) = 3/(X 2 (t)(A - k 2 A b )) + exp(-|z) (d x {t) + d 2 {t)z) . 

Any of these solutions will satisfy the first four Einstein's equations. In order that the rest of the 
Einstein's equations are also valid in the bulk, we must generally have non-trivial components of 
the energy momentum tensor i.e. T\ = T}(t, z), i = 4, D — 1. 



3.2 Bowl-models 

Of particular interest are models with only a single brane with matter, hence one does not have to 
be concerned with the properties and implications of matter on a hidden brane. This means that 
either we have no second brane at all or it only has tension and no matter. Clearly, if we wish to 
have only a single brane, we must have more than one warped extra dimension so that we can wrap 
the brane along a compact extra dimension. The coordinate point that we have to be concerned 
with in such a construction is the origin, z = 0. Having a vanishing stress tensor at z = means 
that the jumps of R and hi must vanish at the origin i.e. the first derivative of the functions must 
be continuous. Clearly, the exponential Ansdtze for the 6j:s do not satisfy these conditions and in 
general it is not easy to see whether solutions where R'(t,0) = (i , 0) = 0, exist or not. 

As a special case we can construct a solution by choosing an essentially common Ansatz for the 
compact extra dimensions: 

bi(z) = bi(0) cosher (kz). (12) 
With this Ansatz, we can solve for R(z,t): 

R 2 (t,z) 



f3\ 



2 



+d 1 (t)e^ 2 ^ a ^ 2 F 1 ( v — -,^,1- V o; „ ,e 2kz ) (13) 



+d 2 (t)e { 



V" 2 - 


4/3 


2k 




4/3 


2k 



2k '2/fc' 

where a = Jjrj-, P = — k 2 A b ) and 2 F\(a, b; c; z) is the hypergeometric function. One can now 

choose one of the (ij(t)'s in such a way that R'(0) vanishes. 

Unfortunately, this solution is not very intuitive, and so, in order to study the properties of the 
cosmological bowl solutions, we restrict ourselves to the case k 2 = 2^^k 2 A b . In six dimensions 
(d = 1) this implies 

R 2 (t, z) = x2 ^ k2Ab Hcosh(kz)) + d(t), (14) 



whereas for d > 1 we write 



R 2 (t,z) = d^t) cosh" (kz) - A 2^ fc2t/ ' ( 15 ) 
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where the exponent is given by v = 1 — ^n-. Note, that the condition R'(0) = has already been 
exploited so that there is only one brane in the scenario which, technically, removes one of arbitrary 
functions d±, (L2 of the Eq. (jl.Sft . Note also, that these solutions correspond necessarily to a positive 
bulk cosmological constant. This model will be studied in more detail in Section 



3.3 Jump equations 

Next we move forward to study the properties of the models on the brane(s). We relate the 
energy momentum tensor of the D-l -dimensional brane to the solutions of D-dimensional Einstein 
equations. We idealize the brane to be a zero width object in z-direction regardless how the brane 
has been formed. Thus, any brane locates at an incontinuity of the z-derivative of the metric tensor, 
in particular at the edge of the space. 

So, the energy content of the brane is related to the global solutions through the jump equations 
The jump equations in the D-dimensional case can be read from the Einstein's equations as 
usual. The brane contribution to the energy-momentum tensor is supposed to be confined to some 
z = z', i.e. it is assumed to have the form 

(Tf ) Br = S(z - z') Tf = 5(z - z') diag(A + /0 , A-p, A-p, A-p, 0, Ai -Pi, A D _ 5 - Pd- 5 )- (16) 

Thus there is a uniform matter distribution also in the directions of the compact dimensions. 
Solving for the different components we get (for general 6j = bi(t, z)): 



R 

M 
77 



D - 

[R!_ 
R 



(D 



D-5 

6)(A + p) + ^ (Ai- Pi 

i=l 



((D-5)(p + p) 
K 2 {p + p) 

^L-(-3(p + p)-(D-6)(A + p) + J2(^-Pi 



(17) 



D-5 



i=l 



D-5 



2(D 



■(4A + p-3p+ ]T (A i -P i ) + (3-L>)(A,--P j 



From these, we see that there are two special cases, D = 5 where there is no pressure (p) (nor 
compact dimensions) dependence on the jump of R, and D = 6, where the brane tension A does 
not contribute either to the jump of R or 77. 



One of the continuity equations (T 0A . A = 0), takes on the brane the usual form 

p + 3H(p + p) = 0, 



(18) 



where H = R/ 'R\ z >. Because we have the freedom to choose A, we may use the standard time at a 
brane by taking X(t) = 1/R(t,z'). Note, however, that if there is more than one brane, standard 
time can be chosen only at one of these. The other constraints can be obtained by noting that 
from the solution of rj, it follows that tj'/tj = R' /R + R/Rdt(R' /R). On the brane, this relates the 
jumps by 



77 



R 



+ H 8 < 



Ml 



[Rf] 
R 



(19) 



=0 



Using the jump conditions and the continuity equation l|18j). we get a relation that binds the brane 
energy density, p, and brane pressures, p, Pi by 



(D-5)(p--p) 



D-5 
i=l 



0. 



(20) 
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Hence, if the jump condition for R' is satisfied and Eq. (|20j) holds, the jump condition for r{ is 
automatically satisfied. Note that in cosmological 5D-models 011)], relation ()19|) reduces to the 
usual continuity equation and hence a solution which satisfies the jump condition for R' will also 
satisfy the condition for jf , which was noted in It is clear that in higher dimensional models, 
some energy can be flowing in/out of the compact dimensions when the universe is not radiation 
dominated (p ^ p/3). 

By assuming that the compact dimensions are static, bj = 0, we get by differentiating the jump 
equation for b'j/bj with respect to time, an equation for each j = 1, d that relates the energy 
and pressures in different dimensions on the brane, 

D-5 

P-3P- E Pi ~ (3 - D)Pj = 0. (21) 

*=i> #3 

Note that these equations and Eq. (|2T)|) are not linearly independent, so that the system is not 
over constrained. By summing the equations for each j, we get again Eq. (|2Uj) and hence we only 
have d equations for d unknowns. 



3.4 Friedmann equations 

With the global bulk solutions and the jump conditions, we can now calculate the Friedmann 
equation on the brane in each case. At first we shall be concerned with models with two branes 
separated by distance L, which may be negative. The branes are located at z = z$ and z = Z\ = 
zq + L. We normalize the scale factor R(t, z) so that R(t, zq) = ao(t) i.e. we live on the brane at 
zq. The Friedmann equation in each case can be calculated by first calculating the jump of the first 
derivative at the two branes, which allows us to express \{t) in terms of cij = — [R']/R\ Zi :s. Since 
the Hubble constant on our brane is given by H 2 = (a'o/ao) 2 = l/(agA 2 ), we can then express H 
as a function of <7j:s. 

The Friedmann equations for the four types of the exponential models are: 



A) H 2 



A + k 2 A b B 2 + AB (<7q + <ri) - 4 ( 7 2 + 2 ^ +g 7 (<r - a x ) - 4q a x 

3 B 2 + ABo l -A 1 (7 -2 ((e 2 ^ L + 1) -2e LE± ? 1± ^ a^e 2 ^ - l)" 1 

tj\ H 2 2(e BL -lVoa 1 -i?(ao-a 1 e BL ) 
] 2(e BL - 1)<ti + B{e BL - 1 - 1o x Le BL ) { ' 

C)H - L(l-a lL ) 

2 B 2 B 2 L + 4cji(2 + BL) + 4a (BL - 2 + Aa x L) 



D) H 



B 2 L + 4cji(2 + BL- 2e BL l 2 ) 



where 7 = \/9B 2 — 24 (A — k 2 Kb). By using the jump equations, one can write these in terms of 
the quantities defined on the branes. 

For the bowl-model the Friedmann equation has a remarkably simple form. We obtain 

k 

H 2 = -[2a coth(/cz ) - ku], (23) 

which is valid for all d > 1. 

We study these equations in more detail in the next Section, where the six-dimensional case is 
considered. Generalization to higher dimensions can then be done straightforwardly. 
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4 Six-dimensional case 



As we have seen from the general expressions, in six dimensions there is no contribution to the jumps 
of R or rj from the brane tension. Also, with only one compact dimension, the extra dimension can 
be identified with IR 2 and we can write the cosmological solution for a model with only a single 
brane wrapped around the origin with no additional matter or tension carrying brane. In this 
Section we study these properties in more detail. 

The (OO)-component of the Einstein's equations in six dimension, after substituting r](t, z) = 
\(t)R(i,z), is 

3 b" b 1 R' , R' . 9 R 9 , . . 

+ IT + 3 T ^ + 3(^) 2 + 3^ = k 2 A b . (24) 



X 2 R 2 b b R v R ' R 
Let us first look for solutions with the exponential Ansatz, b(t,z) = b^exp^kz). Again, de- 
pending on the value of the parameters we have different types of non-trivial solutions, which can 
be read out from Eqs. However, the case C is possible only if k = A B = 0, corresponding 

constant metric component b and the solution given in Eq. IjllJ) . Thus, because A = B 2 = k 2 , the 
cases that remain are: 

A) k 2 / k 2 Kb, k 2 ^ 8k 2 A b /5: 



R 2 (t,z) = di(t) exp((-(-3fc + V24k 2 A b - 15A: 2 )) + 
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d 2 (t) exp(-(-3fc - V24k 2 A b - 15k 2 )) + 



B) k 2 = k 2 A b ^ 0: 

D) k 2 = -8k 2 A b /5 / 0: 



6 V v D " (k 2 -K 2 A B )X 2 



R 2 {t,z)=d 1 (t)e- kz + d 2 {t) + ^ (25) 



R2 ^ ,Z) = ~ \ 2 (t)K 2 A B + exp( ~^ z) + d2 ^ 



In addition to these, we also have the bowl-solutions for which b'(0) vanishes. In six dimensions 
the Ansatz for b is b(z) = cosh(A;z) and generally we write 

,kz , „ „,1 1, , 1 



2kz 

x, —f 

,kz, „ „ A 1, , 1 



R\t,z) = d 1 (t)exp(—(l-a)) 2 F 1 (-,-(l-a),l--a,-e 2kz ) + 



d 2 (t) exp(y (1 + a)) 2 F x q, |(1 + a), 1 + -e 2kz ) + A2(fc2 _ ^ , (26) 



where a = y 1 — |(1 — K ' 2 j^ B ) and 2 Fi(a, b; c; z) is the hypergeometric function. In particular, if 
k 2 = k 2 A b we obtain the special solution given by Eq. (J14[). 



4.1 Jumps 



The cosmological solutions obtained by placing branes at z = zq and z = z\ can now be studied as 
previously. From the jump equations (fTTj) we get: 

,.-2 



[V] 
b 

m 

R 

M 

V 



-(4A i + p i -3p i -3A^ + 3P i ; 



-{pi+Pi + Al-Pi] 



— (Spi-Spi + Af-Pi), 



(27) 
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where Aj and Af denote the tensions of usual four dimensional Minkowski manifold and 9 direction, 
correspondingly. Note how the brane tension, Aj, does not contribute to the jump of R'. The 
exponential Ansatz obviously relates the energy density on the two branes since [6']/6]^ = ±k = 
— [&']/& | ZD where the sign in front of k is same as the sign of the difference Z\ — zq. 

An interesting case may be found if the other brane is assumed to be empty of matter, i.e. 
Pi = Pi = Pi = and even A® = 0. In this case there is only a four-dimensional cosmological 
constant at z\ and o\ = 0. Also we find that k = k 2 A z= q and , in D -case, the negative cosmological 
constant of the bulk reads Ab = — |k 2 A 2 . 



4.2 Friedmann equations 

The Friedmann equations for the cases A and B and D can be read from the previously obtained 
solutions (|22|). again with A = k 2 , B = k so, that 7 = \/2Ak' 2 Ab — 15k 2 . These special cases of 
the equations are rather lengthy and not very intuitive and hence they are not explicitly presented 
here. 

The Friedmann equations on our brane at z = zq have interesting limits when the two branes 
are well separated from each other. In case B we see that 

Lk^>l,H 2 « ^ po+po + A e -P )- A 



K 2 



-Lk%>l,H 2 « \k\—(p + Po + A ti -P ) (28) 

where we have assumed that |fc| ~ |(7j|. Hence, in these limit cases the Friedmann equation is linear 
in pq and is not coupled to the matter on the other brane. Note also that Eqs. 1)22(1 are linear in a$ 
when u\ = recreating the conventional evolution of the scale factor. In the bowl-model, we again 
follow the same procedure as before and find the Friedmann equation on the brane (note, that now 
we only have a single brane): 

K 3 

H 2 = — ^XE coth{kz ){p +P0 + A 9 Q -P ). (29) 

Finally we wish to give the relations using the observed 4-dimensional quantities. By using Eq. 
(|20[) we can express the pressure Po in terms of the 5-dimensional energy density: By integration we 
obtain pi = 3(pi — Pi) + C, where C is a constant. Moreover, in order to recover the usual Friedmann 
equation, we need to express the energy density in terms of the four-dimensional quantity. Because 
the components of the 5-dimensional energy-momentum tensor are related to the 4-dimensional 
counterpart by (T^)4D = J dz d6b{z)5(z — zq) T^. We obtain for the exponential models {Tj[)±d = 
2vr6(0)e fc20 ff and for the bowl model (Tf ) iD = 2vr6(0) cosh(fcz )rf . 

Thus the Friedman equations (|28|) and (|29|) may be written in terms of 4-dimensional quantities. 
We find that they are of the form 

H 2 = ^^(P0 + A eff ), (30) 

where the effective mass scale M e ff and effective cosmological constant A e ff depend on the model. 
In particular, for the bowl model 

8VT .3 ^ 



SM 2 ff 2 2vr6(0) sinh(fczo) 



and A e ff = A 9 — |C It is, however, clear, that the effective scale M e jf is a subject of constraints 
appearing from standard cosmology, nucleosynthesis, recombination etc. because Hubble rate may 
change. Therefore M e ff has to be close enough to Mpi. 
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As the eq. JBHJ gives a relation between the parameter of the model and observed quantities, 
so does the equation, which determines the strength of gravitational interaction. Namely, the 
4-dimensional gravitational constant G4 = Mpf is now given by 

—Rfat) -Q-Jdz , (32) 

where Gq = k 2 /8tt = M~ A . The requirement is thus that the Planck mass Mpi equals now the 
usual value observed in the gravitational experiments. 



5 Conclusions 

In this paper we have studied the cosmological evolution of brane world models with a warped and 
a number of compact dimensions. By studying how the bulk evolves with time we find that the 
induced evolution on the brane(s) is quite naturally linearly proportional to the energy density on 
the brane. Moreover, in a number cases, the usual Friedmann equation is reproduced without any 
dependence on the energy density on the so called hidden brane. 

A particularly nice scenario is the bowl model where only one brane is present and hence there 
is no hidden brane contribution. The standard Friedmann equation is recreated with an effective 
cosmological constant on the brane. Interesting models can also be realized by assuming that the 
other brane carries only brane tension along the non-compact dimensions since then the Friedmann 
equation again has the usual form. Furthermore, standard cosmological evolution can be reached 
in some two-brane models by letting the two branes be well separated from each other. 

An obvious omission in the analysis is the question of radion stabilization. Here we have to, 
in addition to the stabilization of the distance between the brane, be concerned with stabilization 
of the extra compact dimensions. Such considerations are clearly needed when one wishes to fit 
the models considered here into a larger theoretical framework. Here our main concern has been 
to consider these models from a purely phenomenological point of view and see how cosmological 
evolution is affected. In this paper we have also omitted the analysis of the graviton spectrum in 
each case which would be needed in order to study the correction to the Newtonian potential. In 
the six-dimensional case this has been done in 1 1 3 j . 

With both warped and, possibly large, compact extra dimensions, one also needs to keep a 
number of experimental constraints in mind. The appearance of the KK-tower of standard model 
particles sets an upper limit for the size of the extra compact dimensions on the brane, 6j(0) _1 > 
1 TeV. The corrections to the Newtonian potential must also be small enough. In order to alleviate 
the hierarchy problem, the size of the extra dimensional volume must also be large enough, so that 
the effective gravitational constant on our brane is small as given by Eq. 1)32(1 . In six dimensions 
these experimental constraints were analyzed in detail in JS] and these results are probably also 
essentially valid for more than six dimensions. Anyway, we have found that, at this level, adding 
compact dimensions surely increases the number of viable models, despite the new experimental 
constraints set by the KK-particle spectrum. 

Brane world models with warped and compact dimensions are an interesting possibility whose 
implications for e.g. particle phenomenology are not yet fully known. Cosmologically they seem 
appealing since the standard Friedmann evolution with linear dependence on the energy density 
and no appearance of hidden brane quantities, is reached in many cases naturally without any fine 
tuning. 
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